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Abstract 

The occurrence of the El Nino Southern Oscillation (ENSO) has been studied from numerous air-ocean interaction aspects, 
yielding theories implying a positive-only correlation between the time rate of change of sea surface temperature anomaly (0T's/ot) and 
the corresponding change in the depth of the ocean mixed layer (h’). However, an alternate proposal by Chu (1991a,b, 1992) holds 
that there are both positive and negative correlations between dT's/dt and h’. Chu's proposal, based upon the modified Kraus-Tumer 
ocean mixed layer (OML) thermodynamic model, goes further to say that the regions occupied by these positive and negative 
correlation values correspond roughly to those occupied by positive and negative values of a surface forcing function, P, that is an 
indication of the strength of the counteracting surface forcings of wind work and buoyant damping. The model-generated fields of 
surface wind stress, net heat flux, sea surface temperature, and primary (mixed) layer depth are analyzed for a two year ENSO event, 
1986 - 87, covering the region of the equatorial Pacific, 10N to 10S and to 140E to 80W, in order to determine the validity of these 
theories. The results shows a fairly uniform P-field over the entire period, consisting of a negative P-field (indicative of a shallowing 
regime) in most of the western Pacific and in the far eastern part of the eastern Pacific with a positive P-field (indicative of an 
entrainment regime) in between. The correlations between dT's/dt and h' show that there exists both positive and negative correlations 
between OT's/at and h', however, there is no appreciable similarity between the regions occupied by positive and negative correlations 
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ABSTRACT 


The occurrence of the E] Nino Southern Oscillation (ENSO) has been studied 
from numerous aif-ocean interaction aspects, yielding theories implying a 
positive-only correlation between the time rate of change of sea surface 
temperature anomaly (0T's/ot) and the corresponding change in the depth of the 
ocean mixed layer (h'). However, an alternate proposal by Chu (1991a,b, 1992) 
holds that there are both positive and negative correlations between OT's/dt and 
h’. Chu's proposal, based upon the modified Kraus-Turner ocean mixed layer 
(OML) thermodynamic model, goes further to say that the regions occupied by 
these positive and negative correlation values correspond roughly to those 
occupied by positive and negative values of a surface forcing function, P, that is 
an indication of the strength of the counteracting surface forcings of wind work 
and buoyant damping. The model-generated fields of surface wind stress, net 
heat flux, sea surface temperature, and primary (mixed) layer depth are analyzed 
for a two year ENSO event, 1986 - 87, covering the region of the equatorial 
Pacific, 1ON to 10S and to 140E to 80W, in order to determine the validity of 
these theories. The results shows a fairly uniform P-field over the entire period, 
consisting of a negative P-field (indicative of a shallowing regime) in most of the 
western Pacific and in the far eastern part of the eastern Pacific with a positive P- 
field (indicative of an entrainment regime) in between. The correlations between 
OT's/dt and h' show that there exists both positive and negative correlations 
between OT's/dt and h’, however, there is no appreciable similarity between the 
regions occupied by positive and negative correlations and those occupied by 


positive and negative P-values. 
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I. INTRODUCTION 


A. RATIONALE FOR STUDY 

Currently, there are several significant theories containing simple ocean 
thermodynamics that have been developed to study the phenomenon of the El 
Nino Southern Oscillation (ENSO). Certain theories differ, however, concerning 
the correlation between the time rate of change of sea surface temperature 
anomaly (0T';/dt) and the anomaly of upper layer thickness (1), or mixed layer 
depth anomaly (h'). The relationship between 1) and h' will be explained in the 
following background section, but for now, it is sufficient to say that they are 
interchangeable. 

Numerous published theories (Wyrtki, 1975; McCreary, 1983; Cane and 
Zebiak, 1985; Battisti, 1988,1989; and Hirst, 1986 among others) show a positive- 
only correlation between OT',/dt and ). The modified Kraus-Turner (K-T) ocean 
mixed layer (OML) thermodynamic model proposed by Chu (1991a,b, 1992) 
holds, conversely, that there are both positive and negative correlations between 
these two parameters depending on the physical process being considered. The 
purpose of this study is to analyze data to establish the proper correlation 
between these parameters and thus establish a correct basis for the study of the El 
Nino phenomenon. 

1. Upwelling Physics 

The theories based on upwelling physics are repesented by the follwing 
thermodynamic equation (Hirst, 1986): 
oT, /ot+ueTx=k(on-T',),k>0,0>0 (1) 


which implies a positive correlation between OT’;/dt and n’. In the limit of large k, 
this equation reduces to TI’; = On, indicating a Bein veronlt correlation between 
T's and 1). 

2. Ocean Mixed Layer Physics 

Chu (1991a,b, 1992) draws a distinction between “upwelling physics" 
and "OML physics" in dealing with the correlation between OT',/dt and h’. When 
considering "upwelling physics," based upon a cooling upper layer where the 
outgoing heat flux from the base of the mixed layer is greater than the incoming 
heat flux, the following relationship is valid: 

OT',, 0t~h',, (2) 
where there is a positive correlation between QT's / dt and h’. Based upon this 
type of physics, equation (1) of the Hirst model is valid. However, when 
considering "OML physics", where the upper layer may be under either strong or 
weak wind forcing, there exist both positive and negative correlations between 
these two parameters 

OT, /ot~(4+/-)h', (3) 
where there exists a positive correlation for strong wind forcing and negative 
correlation for weak wind forcing. Chu takes this relationship one step further, 
relating these positive and negative correlations to corresponding positive and 
negative values of a surface forcing function, P. This function will be discussed 


further in the following section. 


B. BACKGROUND 
1. Relationship of Mixed Layer Depth to Upper Layer Thickness 
Anomalies 
The value of the upper layer thickness anomaly, 1, is based on the 
position of the main thermocline, and is defined as: 
ii ead (4) 
where g' = gAp / p, 1) is the instantaneous change in the depth of the thermocline, 
and d is the initial depth of the mixed layer. Figure 1 illustrates the physical 


relationship between h and 0. 


Convergence Go 


d Ux 





NN 


| 


Figure 1. Illustration of upper layer anomaly, 1), where d represents 
the initial thickness, and h represents the instantaneous mixed layer 
depth (from Yamagata, 1985). 


Essentially, N represents the change in the depth of the thermocline 
lying just below the mixed layer, and therefore corresponds to the mixed layer 
depth anomaly, h’. 

2. Significance of Surface Forcing function, P 
A discussion of P should begin with the modified K-T model that 


introduced P as a surface forcing function. The value of P, whether positive or 


negative, is indicative of the predominant surface forcing, wind work or buoyant 


damping respectively, on the OML, and is defined by the following equation: 
la Cu. ‘a C2 Agh,Qp ! PwC pw (5) 


(a) (b) 

where C, and C> are tuning coefficients, uy» is the water surface friction velocity 
defined as (CpPao/Pwo)!/2 | Vp], Vp is the surface wind velocity, hy is the 
thickness of the OML, Qo is the net surface heat flux, pag and Pwo are 
characteristic values of surface air density and water density, and @ is the thermal 
expansion coefficient. The values used for the tuning coefficients, C; and Co, are 
1.0 and 0.5 respecitively (Schneider and Muller, 1990). Term (a) represents the 
surface forcing due to wind work, and term (b) defines the surface forcing due to 
buoyant forcing/damping. Wind work acts to increase turbulence in the mixed 
layer, and thus increase the depth of the mixed layer. Buoyant damping is caused 
by downward heat flux, +Qo, which acts to stabilize the mixed layer by heating 
the upper portion of the mixed layer leading to stratification and thus a 
shallowing of the mixed layer. Strong surface wind forcing, uy*, will result in a 
positive value for P. Likewise, strong downward heat flux, Qo, accompanied by 
moderate to weak wind forcing, will result in a negative value for P. Generally, 
the value of P varies from positive to negative throughout the ocean basin based 
on the diurnal and seasonal variations in Qo and the synoptic and seasonal 
variations in uw*. 


The modified K-T model incorporates this P-value in the following 


equation: 
aT", Jatt tiy eVT, = Aw PIS a (6) 
Pwol pwhy 
(a) (b) (c) 


where term (a) is the rate of change for Ts, term (b) the horizontal advection of Ts, 
and term (c) the thermodynamic forcing, including the heat flux at the base of the 
mixed layer, Q-p. Of significance in term (c) is the heavyside step function, A\(P). 
This function is referred to as an OML "switcher." For positive values of P, the 
value of the switcher is 1, and for zero or negative values of P, the value of the 
switcher is 0. These two cases have the following effect on the modified K-T 


model in equation (6): 


a. P<O 
In this case, the value of A(P) is 0, yielding the following equation: 
aT, /at+ii, «VT, =—<2— (7) 
Pw0°pwhy 


Here, a negative P-value results from weak wind forcing and/or strong downward 
heat flux, causing a shallowing in the OML. Since the mixed layer will heat more 
rapidly as its depth decreases, there will be a positive correlation between h' and 


OT';ot. This case is illustrated in figure 2. 
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Figure 2. Kraus-Turner OML Thermodynamics under weak wind 
forcing, P<0 (from Chu, 1991a). 


b P>O 
In this case, the value of A(P) is 1, yielding: 
oT’, /ot+u,, ° VT, -Zo~ Gn (8) 


i Oe pwilw 


As noted before, a positive P-value results from strong wind forcing, which causes 
entrainment of deep water into the OML, as indicated by the presence of the Q.y 
term. Since Q.n is greater than Qo, the OML is referred to as a “cooling layer” 
(Miller, 1976), and an increase in the mixed layer depth will be accompanied by a 
corresponding increase in OT';/dt. Thus, there would be a positive correlation 
between h' and OT',/ot. In physical terms, the thicker the layer becomes the less 


rapidly it will cool. This case is illustrated in figure 3. 
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Figure 3. Kraus-Turner OML thermodynamics under strong wind 
forcing, P>0 (from Chu, 1991a). 


Il. DATA COLLECTION AND PROCESSING 


A. DATA COLLECTION 
In order to effectively study the correlation between dT',/dt and h' and its 
relationship to the P-value, the following data sets are required: monthly mean 
net surface heat flux (Qnez), monthly mean surface wind stress (Ty), daily mean sea 
surface temperature(T,), and daily mean mixed layer depth (h). The area of 
coverage includes the majority of the equatorial Pacific. Therefore, the study area 
is chosen as 140 E to 80 W and ION to 10S inclusive. Figure 4 is a map of the 
equatorial Pacific region. The time period is chosen to span a designated ENSO 
event, January 1986 - December 1987. Data was collected for this study from 
two sources, Fleet Numerical Oceanography Center (FNOC) in Monterey, 
California, and the National Center for Atmospheric Research (NCAR) in Boulder, 
Colorado. The data products received from FNOC, T; and h, were generated by 
the Thermodynamic Ocean Prediction System (TOPS) - coupled Expanded Ocean 
Thermal Structure (EOTS) Analysis model. The data products received from 
NCAR were Qnet, Tw, and Ts, and were archived from the European Center for 
Medium Range Weather Forecasting (ECMWEF) synoptic forecasts. A 
description of these data sets and their forecasting and modeling methods is 
provided below. 
1. FNOC Archive Data Sets 
Data was acquired through FNOC from the archives of the Data Base 
Management Division. During the time period of 1986-87, the data products of 
sea surface temperature and primary (mixed) layer depth were generated by the 


TOPS- coupled EOTS (TEOTS) analysis technique. A discussion of the FNOC 


data will include a section on TEOTS analysis methods and another section on 


TEOTS data product output. 





Figure 4. Equatorial Pacific region 
a. TEOTS Analysis Method 

The EOTS analysis scheme combines ocean thermal climatology 
with real time thermal observations, applying a series of corrections via the Fields- 
by-Information-Blending (FIB) analysis technique. The ocean climatology is 
drawn from an archive of coarse resolution (320 km) monthly climatology. The 
observations are provided by real time ship and satellite data. Daily Ship 
observations consist of about 150-200 XBT's , and 1200-2000 ship injection/ 
bucket temperature measurements. Daily Satellite measurements consist of 
50,000-80,000 Multi-Channel Sea Surface Temperature (MCSST) retrievals from 
NOAA. This system is limited, however, since local bursts of atmospheric forcing, 
which cause considerable change in the mixed layer depth and sea surface 
temperature are often neglected by using such a sparse number of real time 
measurements. For this reason, FNOC adapted a system, endorsed by Elsberry 


and Garwood (1980) and Clancy and Martin (1981), coupling a synoptic mixed 


layer model (TOPS) to the basic EOTS analysis to generate a continuously 
updated, real-time, analysis-forecast-analysis system. The gist of the system is to 
initialize the TOPS model daily using the EOTS analysis fields, produce a forecast 
field driven by observed and forecasted atmospheric forcing, and subsequently 
feed this field back into the EOTS analysis as a first-guess for the following 
analysis day. This system is collectively referred to as TEOTS, and appears in the 
form of a flow chart in figure 5 (FNOC, 1987). 
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Figure 5. TEOTS Flow Chart (from FNOC, 1987) 

One concern in deciding whether or not to use the TEOTS- 
generated data involves the relationship between Ts and h. One of main goals of 
this study is to determine a correlation between these two parameters, and 


therefore it is necessary to establish that there are no physical assumptions 


concerning this correlation within the TOPS forecast that might bias the results of 
this study. The following is a brief discussion of the physics used in TOPS. 

The TOPS consists of four prognostic equations representing the 
conservation of temperature, salinity, and momentum in the upper ocean (FNOC, 


1987). 
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(10) 
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where T is the temperature, S is the salinity, u, v, and w are the x, y, and z 
components of the current velocity, F is the solar radiation flux (positive 
downward), Py is the reference density, c is the specific heat of seawater, D is a 
damping coefficient, v is a diffusion coefficient, z is the vertical coordinate 
(positive upward from the surface), and ug, Va, and wa are the x, y, and z 


components of the "advection" current. The parameter A is a horizontal eddy 
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diffusion coefficient whose value depends on the horizontal grid spacing. 
Equations (9) - (12) exclude the effects of geostrophic flow, which was 
determined to have a negligible effect on mixed layer dynamics. 

These equations were parameterized by the following boundary 


conditions (Mellor and Durbin, 1975 and Clancy and Martin, 1981): 


w'T' =—qS; a (13) 
w' S' =~LqSy, 8 (14) 
wiu' =—lqSy a (15) 
wv =—LgSy a (16) 


where @ is the turbulence length scale, q is the square of twice the turbulent 


kinetic energy and is defined by the following equation: 


au (avy g do q° 
eqS — — éqS1,| —-—— |-—- = 0 17 
‘ (2) (=) } : »( 2 | 1s mn 


The three terms in equation (17) represent wind forcing, or shear production, 
buoyant damping, and viscous dissipation respectively. The variables Syy and Sy 


are functions of the gradient Richardson number, Ri, defined as follows: 
_& 9p 


Re Pw OZ (18) 


-@| 


il 


The Richardson number basically indicates the predominant forcing between 
Shear production and buoyant damping. If the Richardson number exceeds a 
critical number, Ri = .23, buoyant damping stratifies the mixed layer sufficiently 
enough to inhibit the turbulent mixing of shear production, and the layer enters a 
shallowing regime (Miles, 1961 and LeBlond and Mysak, 1978). Conversely, if Ri 
< .23, shear production dominates, resulting in entrainment and a deepening of 


the mixed layer. The extent to which the mixed layer deepens is dependent on 


the magnitude of the turbulent length scale, 2, defined as follows: 
0 
Ie | z\gdz 
0 
ie qdz 


As seen in the progression through the prognostic equations (9) - (12) to 


£=0.1 (19) 


the closure of equation (19), there are no assumptions made regarding the direct 
relationship of mixed layer depth to sea surface temperature. The TOPS model is a 
very meticulous mixed layer forecasting system whose only assumption excludes 
the presence of geostrophic flow. For this reason, the TEOTS analysis system was 
determined to be unbiased regarding the objective of this study. 
b. TEOTS Analysis Output Description 

The output of the TEOTS analysis took the form of two different 
grid schemes, one for T,; and one for h. Mixed layer depth is output on a grid 
using a polar stereographic projection. This type of scheme consists of 63 X 63 
grids overlaid on maps of the Northern and Southern Hemispheres whose origins 
are the poles. The grid spacing is approximately 320 kilometers. A disadvantage 
of this grid scheme, when dealing with equatorial data, is that the resolution of the 
data decreases radially outward from the polar origin, yielding the worst 


resolution at the equator. Sea surface temperature is output on a grid using a 


he 


spherical, or Mercator, plot. This type of scheme consists of a 73 X 144 grid 
overlaid on a Mercator projection of the globe. The grid spacing for the spherical 
plot is approximately 2.5°, or 280 kilometers. The contoured plots of monthly 
mean mixed layer depth are found in Appendix D. The contoured plots of 
monthly sea surface temperature are found in Appendix C. 
2. NCAR/ECMWE Archived Data Sets 
As previously mentioned, the ECMWF data sets were received from 
NCAR. A description of the numerical models used to generate this data will be 
quite brief due to the emphasis on mixed layer dynamics rather than the physics 
of the processes driving atmospheric forcing. 
a. ECMWF Numerical Weather Prediction 
The wind and surface flux fields provided by ECMWE are 
generated by an operational model using numerical weather prediction. The 
model used for wind analysis and prediction produces 6-hourly forecasts on 15 
standard levels in the vertical and on a 1.125° horizontal grid. The model 
combines observations obtained within a six hour window of the analysis time 
with a multivariate height, wind, and thickness system (Lorenc, 1981) from the 
preceding analysis period to produce a forecast for the following period. The 
method used to combine these two inputs, Statistical interpolation, assigns 
weights to the these inputs, and uses geostrophy and other error checking criteria 
to smooth the results (Lonnberg, 1982). A system flow chart is shown in figure 8. 
Once the surface wind field has been generated, the surface wind 


stress fields are computed using the relation: 
T =pCpU 2 (20) 
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where p is the density of air in the planetary boundary layer, Cp is the drag 
coefficient over the surface, and U is the 10 meter wind provided by the 


foregoing model. 
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Figure 6. ECMWF analysis / forecast flowchart 


Surface flux fields are simply a byproduct of the standard outputs (winds, 
temperature, potential temperature, and humidity) of the ECMWF model. 
Shortwave and longwave radiation are computed using straightforward 
atmospheric physics combined with satellite cloud imagery to form an algorithm 
modeling the behavior of the radiation budget. Radiation affects the ocean 
mixed layer both directly, by increasing or decreasing the heat flux driving the 
buoyant damping on the layer, and indirectly, by affecting the development of 
clouds and weather patterns which have a local effect on shear production. A 


schematic of these radiative processes is shown in figure 9. Longwave radiation 
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emitted by the surface is computed using the modified Stefan-Boltzman law 

R, =—-e0T ‘ (21) 
where € is the emissivity of water, 6 is the Stefan-Boltzman constant, and T is the 
temperature of the surface in degrees K. Since longwave radiation is emitted not 
only by the surface but by clouds and the atmosphere, the estimation of 
longwave radiation is a fairly complicated function of the earth and satellite- 
based soundings. Shortwave radiation at the top of the atmosphere 1s computed 
using an inverse square law coupled with Lambert's cosine law to yield 

Ro = So (Gn /d)’ cosy (22) 

where So is the solar constant, d,,/d is the ratio of the mean to the 
actual distance from the sun, and y is the solar zenith angle. In order to obtain the 
actual surface value of shortwave radiation, absorption and scattering 


predominantly by clouds are also considered. 
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Figure 7. Radiative Processes in the Planetary Boundary Layer (from 
Wergen, 1982) 
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Latent heat flux is computed from the basic equation 


H, = LE (23) 
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where Le is the latent heat of evaporation/condensation and E 1s the rate of 
evaporation/condensation determined by the soundings. Sensible heat flux is 
calculated using the bulk aerodynamic approach based on the Monin-Obukhov 
similarity theory. This theory assumes that near the surface, vertical wind and 
potential temperature gradients depend only upon height, surface sensible heat 
flux, surface momentum flux, and the thermal expansion coefficient of air. 
Accordingly, ECMWE uses the temperature and winds at their lowest sounding 
level (33 m) and the surface temperature to compute the sensible heat flux at the 


surface as 


k 4 
2 in(h/Z). ee y a nh 


where k is von Karman's constant, h/zo is the ratio of height to roughness length, 
Vn is the wind velocity at the 33 m level, A@p is the change in potential 
temperature between the surface and the 33 m level, and G(h/zo,Ri) 1s an 
analytical function that is continuously updated based on observations and 
performance of the model (Wergen, 1982). The Richardson number, Ri, is given 


by 
ghd, 


alV, |" 





(25) 


= 


and is an indication of stability within the surface boundary layer. 
b. ECMWEF Output Description 
The ECMWF data fields accessed through NCAR were archived in 
grid-in-binary (GRIB) format. In order to compute the surface friction velocity, 
the following data sets were extracted: u-component of surface wind stress, and 
v-component of surface wind stress. In order to compute the net surface heat 


flux, the following data sets were extracted: net surface shortwave radiation, net 
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Surface longwave radiation, surface sensible heat flux, and surface latent heat 
flux. The data fields are derived from 6 hour forecasts valid for 0 Z and 12 Z, and 
are plotted on a 1.125° grid. The contoured plots of monthly mean net surface 
heat flux are found in Appendix A. The contoured plots of monthly mean water 


surface friction velocity are found in Appendix B. 


B. DATA PROCESSING 

Data processing consisted of two basic steps. First, wind and heat flux data 
fields are combined with mixed layer data to compute values for the surface 
forcing function, P. Sea surface temperature and mixed layer depth were then 
analyzed to determine the correlation between OT',/dt and h’. 

1. P-Value Computation 

The P-value fields were computed by equation (4). Since the majority of 
the data required for this computation came from the ECMWF forecasts, a 1.125° 
grid was chosen to represent the P-field. One P-field was calculated per month 
from January 1986 - December 1987. The following equation was used to 
compute the value of the net surface heat flux, Qo: 

Op = Qs + Qw + Qs + Q, (26) 
where Qgy is the shortwave radiative flux, Qjy is the longwave radiative flux, Qs is 
the sensible heat flux, and Q) is the latent heat flux. All fluxes are in the units of 
watts per meter squared (W/m2). These terms were defined as positive 
downwards and vice versa as archived, and thus the terms were simply added 
together to yield Qo. 

In order to obtain the value of water surface friction velocity, uw* the u 
and v-components of surface wind stress were first combined into a scalar value 


of surface wind stress, in units of N/m2. This value of wind stress was then 
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converted to wind speed in the planetary boundary layer, Vp, in units of meters 
per second (m/s). The following equation was then used to calculate uy*: 
u_ + =(CpPa0/ Pwo)” |VsI (27) 

The only remaining variable in the calculation of P is h. The only 
modification required here was to convert the depth from units of feet to meters. 
However, since the grid spacing for mixed layer depth was irregular and greater 
than that of winds and heat flux, a bivariable interpolation routine was used to 
interpolate the polar stereographic mixed layer depth grid field to the 1.125 
Mercator grid. Contoured plots of monthly mean P-values are shown in 
Appendix E. 

2. Correlation Computation 

The first step in calculating the correlation between OT',/ot and h' was to 
match the grids between T; and h. For ease of plotting, the 2.5° grid for Ts; was 
chosen. Here again a bivarible interpolation routine was used to interpolate the 
polar stereographic mixed layer depth grid field to the spherical Ts grid field. The 
objective of the data processing here is to obtain a correlation field based on daily 
trends of OT';/dt and h' over each of the 24 months. Contoured plots of these 


correlations are shown in Appendix F. 
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III. SURFACE FIELDS 


The daily surface data sets were averaged over each of the 24 months of the 


two year period to yield monthly mean averages. 


A. HEAT BUDGET 
In order to see the time variation of net surface heat flux, we plot the zonal 
values of net surface heat flux versus time for latitudes 5N, 0, and 5S (figure 7). 
1. Western Equatorial Pacific 
The heat budget in the western equatorial Pacific is characterized by a 
downward net surface heat flux throughout 1986-87. The highest downward 
heat flux occurs in the southwestern equatorial Pacific during the period of 
September -November 86 (225 W/m2). In the northwestern equatorial Pacific, the 
downward heat flux values oscillate between strong during the spring and early 
fall to weak during the summer and winter. In the southwestern equatorial 
Pacific, downward heat flux values range between high in the early spring and 
late fall to weak during the summer months. 
2. Mid Equatorial Pacific 
The heat budget in the mid-equatorial Pacific is characterized by a 
predominantly upward net surface heat flux over the course of the two year 
study. The maximum upward heat flux occurs during Aug 86 (255.9 W/m2) in the 
southern portion of the mid equatorial Pacific. North of the equator, heat fluxes 
are strongest during the winter, diminishing elsewhere. South of the equator, 
fluxes are again predominantly strongest in the winter, but proceeding further 


south, fluxes increase strongly in the summer. 
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3. Eastern Equatorial Pacific 
The heat budget in the eastern equatorial Pacific has a typically 
downward net surface heat flux over a period of two years. The maximum 
downward heat flux occurs during Feb 87 (199.2 W/mZ2) in the north. There is, 
however, no distinct seasonal pattern in the temporal fluctuation of heat flux in 


this region. 


B. SURFACE FRICTION VELOCITY 
The water surface friction velocity, uy* as described by equation (5) will be 
referred to as wind speed since it is proportional to wind speed for a constant 
drag coefficient. The time evolution of zonal winds are plotted in the same 
manner as net surface heat flux and are shown in figure 8. Meridional winds are 
not plotted due to their low relative correlation with sea surface temperature 
change as shown by McPhaden and Hayes (1991). 
1. Western Equatorial Pacific 
a. North Equatorial Latitudes 
The western equatorial Pacific is dominated by easterly winds over 
most of the two year period. The maximum easterly wind speed in the north (.26 
cm/s) occur in December 1986. The trend in the north is typically strong easterlies 
during the winter and early spring. 
b. South Equatorial Latitudes 
The maximum easterly winds in the south (.23 cm/s) occur in June 
1986. The trend in the south is typically strong easterlies during the summer 


months diminishing to weak during the winter months. 
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Figure 8. Time evolution of net surface heat flux (Qo) for latitudes (a) 
5N, (b) 0, and (c) 5S across the equatorial Pacific. The value of Qo is 
defined as positive (solid) downwards, and negative (dashed) upwards. 
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2. Mid Equatorial Pacific 
a. North Equatorial Latitudes 
The mid equatorial Pacific is a region of seasonal wind variability. 

In the north the easterly winds dominate in late winter to early spring with the 
maximum wind speed (.36 cm/s) occurring in January 1987, and the westerlies 
dominate in the fall with maximum wind speed (.38 cm/s) occurring in September 
1986. The summer months represent the transition between the easterlies and 
westerlies and thus the near-zero wind speeds are fairly consistent with the 
occurrence of the doldrums. 

b. South Equatorial Latitudes 

In the south, the seasonal shifts between easterlies and westerlies, 
differ from those in the north. Easterlies dominate in the early spring to late 
Summer, with the maximum value (.26 cm/s) occurring in May 1987. The 
transition between easterlies and westerlies occur here in the early fall and late 
spring. The westerlies dominate during the late fall and winter months, with 
maximum value (.28 cm/s) occurring during March 1987. 
3. Eastern Equatorial Pacific 

The eastern equatorial Pacific follows the general patterns and trends of 
the mid equatorial Pacific. In the north, however, the shift between easterlies and 
westerlies occurs earlier, in the spring, and thus the maximum westerlies (.61 cm/s) 
occur in the summer, specifically July 1986. The maximum easterlies (.35 cm/s) 
occur in the early winter, specifically December 1986. In the south, the winds are 


fairly weak over the the two year period. 
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Figure 9. Time evolution of water surface friction velocity (u,*) for 
latitudes (a) 5N, (b) 0, and (c) 5S across the equatorial Pacific. The 
value of uy* is defined as positive (Solid) for westerly winds and negative 
(dashed) for easterly winds. 
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C. SEA SURFACE TEMPERATURE 

The time evolution of sea surface temperature, Ts, is plotted using the time- 
longitude plots used in the previous two sections and is shown in figure 9. 

1. Western Equatorial Pacific 

The western equatorial pacific possesses the highest values of sea 

Surface temperature throughout the equatorial Pacific. These high values 
typically increase to the south. The highest values of sea surface temperature 
occur during the period October-December 1986 and March-July 87 (in excess of 
30° C) in the southern portion of the mid-western equatorial Pacific. Temporally, 
the western equatorial Pacific has by far the lowest sea surface temperature 
variability ( 1° C). These findings are consistent with values obtained from 
Reynolds NMC/CAC Surface Temperature Analysis and the findings of 
McPhaden and Hayes (1991) and McPhaden (1988). 

2. Eastern Equatorial Pacific 

The eastern equatorial Pacific has the lowest values of sea surface 

temperature, generally decreasing to the southeast. The lowest values occur 
during September 1986 and August-September 1987 (about 20° C). The general 
trend of sea surface temperature shows warm values propagating eastward, 
peaking in March-April of each year, and propagating westward, peaking in 
September-October of each year. The temporal variability of temperature 
anomaly was also plotted versus time in the same manner as the latter. The results 
show that the eastern equatorial Pacific has the highest variability in sea surface 
temperature [O(1 C)], with the highest variability occurring in the far southeastern 


equatorial Pacific in March 1986 (3° C) and March 1987 (5° C). 
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Figur olution of sea surface temperature (T;) for 
latitu rhe a iy *(b) 0, and () 5S , across the equatorial Pacific. 
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IV. TWO REGIMES OF THE OML 


The surface forcing function, or P-value is an indication of the dominant 
regime of the ocean mixed layer, whether entrainment (deepening), or shallowing. 
As mentioned previously, the surface forcing impacting these regimes are wind 
shear, or shear production, and heat flux, or buoyancy damping. These fields 
have been described in the chapter III. In addition to these variables, the mixed 


layer depth is also used in the calculation of the P-value. 


A. MIXED LAYER DEPTH 
Mixed layer depth is plotted and analyzed as per the previous sections. The 
plot of time evolution of mixed layer depth versus longitude is shown in figure 
10. 
1. Western Equatorial Pacific 
The western equatorial Pacific has an intermediate range of mixed layer 
values, with higher values increasing to the east. The northern region has the 
largest values of mixed layer depth during the spring of each year (30 m). The 
lowest values of mixed layer depth occur during late summer each year (20 m) 
with values increasing toward the equator. The southern region has the largest 
values of mixed layer depth during late fall of each year (30 m), with the lowest 
values occurring during the spring of each year (17-19 m). 
2. Mid Equatorial Pacific 
The western portion of the mid-equatorial Pacific has largest values of 
mixed layer depth. The highest values of mixed layer depth occur just below the 
equator (50 m) and are sustained throughout the two year period with a 


comparatively low temporal variability (3-4 m). 
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3. Eastern Equatorial Pacific 
The eastern equatorial Pacific maintains the lowest values of mixed layer 
depth, particularly in the far eastern portion. The lowest values (5 m) occur in the 
north (7.5 N) during the spring of each year. The lower latitudes (south of the 
equator) tend to have the highest values of mixed layer depth. The temporal 
trend reveals a deepening mixed layer in the fall and winter, returning to a 
Shallow mixed layer depth during the spring and summer. This trend is more 


pronounced in the region south of the equator. 


B. P-VALUE 
Based on equation (5), the value of the surface forcing function, or P-value, is 
computed using the water surface friction velocity, uy*, the mixed layer depth, h, 
and the net surface heat flux, Qo. A plot of P-value, averaged over two years at 
each grid point, versus time is shown in figure 11 for latitude SN, the equator, and 
D0. 
1. Western Equatorial Pacific 
The western equatorial Pacific is dominated by negative P-values. In 
terms of the mixed layer, this would tend to indicate a shoaling regime. The 
largest negative value (2.3 (cm/s)3) occurs to the far west (150 E) at the equator in 
February 1986. Although there is some temporal variability in the P-value in this 
region, there is no marked trend. 
2. Mid Equatorial Pacific 
The mid equatorial Pacific is characterized by positive P-values. In terms 
of the mixed layer, this would indicate an entrainment regime. The largest 


positive value (2.49 (cm/s)3) occurs in the middle (150 W) just south of the 
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equator in August 1986. There is no temporal trend in the magmitude of the P- 
value. 
3. Eastern Equatorial Pacific 
The eastern equatorial Pacific has both positive and negative P-values. 
The positive values of P occur toward the west, blending in with the positive 
values from the mid equatorial Pacific. The negative values occur to the 
southeast. The largest negative value (-0.58 (cm/s)3) occurs to the far south in 


December 1986. 
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Figure 12. Average P-values for latitudes (a) 1ON through 0, and (b) 0 
through 10S across the equatorial Pacific. 
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V. CORRELATIONS BETWEEN dT's/dT AND H' 


Correlations between OT',/ot and h' are computed with h' leading OT’;/dt by 
one day. Why is this lag time used? McPhaden and Hayes (1991) showed that 
the strongest correlation between sea surface temperature and wind forcing 
occurred with wind forcing leading temperature by one day (figure 11). It can 
therefore be assumed that mixed layer depth, which is driven by turbulence 


created by wind mixing, must therefore lead temperature on the same time scale. 
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Figure 13. Cross correlations between SST and wind speed, zonal wind 
pseudostress, and wind pseudowork. Cross correlation extrema (r) and 
the lag (in days) at which they occur are shown. In each case, the lag 
indicates that cold SST follows high winds by 1 day (from McPhaden and 
Hayes, 1991). 
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In order to study this correlation and its relationship to P-value, the following 
process was used. First, JT';/ot and h' were correlated based on daily values over 
the course of each of the 24 months, yielding monthly correlation fields. These 
fields, shown in Appendix F, are shaded dark for regions having negative 
correlations and light for regions having positive correlations. Comparing these 
fields to the P-fields shown in Appendix E, it is clear on the whole that there is no 
similarity between respective positive and negative regions as would be expected 
by "switching" theory. 

The emphasis now shifts to choosing regions in the equatorial Pacific that 
exhibit typically strong entrainment and shallowing tendencies, based upon 
positive and negative P-values respectively. The regions chosen, based upon 
figure 11, are latitudes (1ON-10S) along 165W for entrainment, and latitiudes 
along 160E for shallowing. The monthly correlations at these points are plotted 
over time, and are shown in figures 12 and 13. The result here is that over time, 
both of these regions undergo positive and negative correlations between OT’,/ot 
and h’, regardless of the region chosen for its entrainment/shallowing tendency. It 
is clear from the apparently random oscillatory behavior of these correlations that 
the correlations between OT’,/ot and h' that we would expect from the modified 
Kraus-Turner model are not dominant on time scales of a day or more. 

It is possible, however, that although the correlations between OT’',/ot and h' 
are not significant on these time scales, that the correlations between simply T's 
and h' might be significant. A test of this hypothesis is meaningful since the 
correlation between these two variables, T', and h' is also commonly accepted as 
being only positive. Correlations between daily values of T'; and h' are computed 


over the two year period at 5S and longitudes 160E and 90W (shoaling regions) 
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and 165W (entrainment region) with the results plotted in figure 16. At 90W, the 
highest correlation (-0.8) occurs with h' leading T’, by several days. The 
difference at this longitude is that instead of a negative maximum correlation for 
h' leading T's, the maximum correlation is positive (.35) and occurs with h' leading 
T's by nearly 100 days. Of importance here is that negative correlations are 
expected between temperature and mixed layer depth in regions of negative P- 
values (90W and 165E) and conversely positive correlations in regions of 


positive P-values (160E). 
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Figure 14. Monthly correlation of OT';/ot and h' versus time for 
longitude 165W and (a) latitudes 10N - 0, (b) latitudes 0-10S. 
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Figure 15. Monthly correlation of 0T',/ot and h' versus time for 
longitude 160E and (a) latitudes LON - 0, (b) latitudes 0-10S. 


35 








ae 
ne 
we 
t 
f f 
p "i 2 eae, ee as Seer 
Q 
' 
a 
/ 
4 
ttn 
uy ti so — ie aie afi = 
emt ray 1m wt a tm temy faut Riu 
fap (tay) 
(hy 
4 ‘ 
a} 
n? 
ad 
5 : 
é t —_— ma —- Fo — im 
9 
at 
1 " 
ih} 
q) 1 Se rao ee a = dom ode oho - 
Re wee sey " pum (0 eam Riw 
fp tefavcd 
(} 
fit e 
net 


ae 


1g oh. ate he -t- - 
wowed Aen mn wwe " Mi fim 


Figure 16. Correlation between daily values of T', and h' for the period 
1986-87 with normalized values of correlation plotted against the 
lag(days) at latitude 5S and longitudes (a) 90W, (b)165W, and (c)160F. 
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VI. SUMMARY AND DISCUSSION 


In studying the correlation between the time rate of change of sea surface 
temperature anomaly, OT';/dt, and the mixed layer depth anomaly, h', several 
findings are apparent. The most significant finding is that there exist both 
positive and negative correlations between QT';/dt and h' at the various grid 
points chosen across the equatorial Pacific. This finding runs contrary to the 
upwelling physics theory showing a positive-only correlation. Since forward 
differencing was used in calculating OT’,/dt, the correlations were based on a one 
day lead of h’. The reason for using this lead time is that mixed layer depth and 
sea surface temperature vary with mesoscale atmospheric variability (Elsberry and 
Garwood, 1978), and mixed layer depth changes caused by wind stress induced 
turbulence lead sea surface temperature changes by about one day (McPhaden 
and Hayes, 1991). 

A second finding is that in the western equatorial Pacific (140E-180E) the 
monthly averaged values of the surface forcing function, P, are typically negative 
throughout the course of the two year period, indicating that this region is 
predominantly a shoaling regime. This is due largely to the strong downward net 
surface heat flux in this region. In the mid-equatorial Pacific (180E-120W), the 
monthly averaged values of the surface forcing function, P, are typically positive 
indicating that this region is predominantly an entrainment regime. Finally, the 
eastern equatorial Pacific is characterized by both positive and negative values, 
with negative values typically to the east indicating a general shallowing regime. 

Although the correlation patterns between OT’,/ot and h' do not fit the P- 


value distribution with any degree of confidence, there is some similarity between 
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the correlations of T’; and h' and P-values. It is found that h' leads T’; by several 
days with a strong negative correlations (-0.8, -0.6) in regions of significantly 
negative P-values (QOW, 160E). The converse is not so conclusive in that h’ leads 
T’, by nearly 100 days with a relatively strong positive correlation (0.35) in a 


region of positive P-values (165W). 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A positive-only correlation between OT';/dt and h' is commonly accepted in 
current thermodynamic theories regarding the El Nino phenomenon. In order for 
these theories to be valid, this relationship must hold true, particularly in the 
equatorial Pacific. Based on the results of this study, this positive-only 
correlation is not true in the equatorial Pacific. It would be premature at this 
point, though, to discount the validity of these theories based solely on the results 
of this study. To more sufficiently study this correlation and its relation to the P- 
value, the following study is recommended. First, the equatorial Pacific should be 
broken up into three regions, western, mid, and eastern, to investigate the unique 
surface forcings occurring in each. In each of these regions, thermistor chains or 
current meters should be utilized at, above (about 5N), and below (about 5S) the 
equator to measure hourly surface and subsurface temperatures. The subsurface 
temperatures could then be analyzed to determine mixed layer depth anomaly 
and subsequently correlated with the time rate of change of sea surface 
temperature anomaly. This finer time scale would not only conform more closely 
to the time scale of turbulence in the mixed layer, but would allow more flexibility 
in adjusting the lead time used in the correlations between mixed layer depth and 
sea surface temperature. The conflicting theories regarding this correlation could 


then be better resolved. 


ay 
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APPENDIX A. MONTHLY NET SURFACE HEAT FLUX 
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APPENDIX B. MONTHLY WATER SURFACE FRICTION 
VELOCITY FIELDS 
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